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SUCCESSIVE APPROXIMATIONS FOR BEAMS 
ELASTIC FOUNDATION 


method resembling the well known moment-area procedure presented 
for the solution beams elastic foundation. exten- 
sion, adaptation, and systematization the famed Vianello-Stodola procedures 
are very briefly mentioned August and Miklés Because 
the inherent power these procedures, reasonably accurate results are ob- 
tained simply and quickly. The difficult mathematics that occur the usual 
analytical solutions are completely avoided. 


INTRODUCTION 


The supporting media, such soil, well load-carrying beam footing 
that rests it, usually display certain amount elastic behavior; yet, for 
simplicity, customary most cases consider the footing infinitely 
rigid. The design based the assumption uniform soil pressure distribu- 
tion the centroid the footing coincides with the resultant the applied 
forces. Similar simplified treatment used for eccentric loads. Such 
assumption becomes less and less accurate the footing the load-carrying 
beam becomes relatively more flexible. For example, the case railroad 
rail supported ties, the ordinary solution would lead erroneous results and 
cannot used. There important class problems which flexibility 
the beam elastic foundation must considered. 

the supporting medium considered true elastic continuum, the 
problem becomes unusually difficult and only few‘ rigorous mathematical 
solutions exist. However, the reaction forces are assumed proportional 
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the deflection the beam that point, the problem becomes greatly simpli- 
fied yet, the solution the simplified problem appears yield results that are 
reasonable agreement with the existing exact ones. the case founda- 
tion problem, rather complex soil phenomena are involved and extreme 
mathematical refinement the solution practical problem seems unwar- 
ranted. Moreover, many problems, such those occurring the stiffening 
rings thin shells, grillage beams, etc., the foregoing assumption holds rigor- 
ously. 

The idealized foundation assumed here may thought one composed 
infinite number independent springs. The elastic constant these 
“springs” represents the foundation modulus. These springs can take both 
tension and compression forces. This assumption was first 
Eduard Winkler. Many subsequent investigators used this assumption, per- 
fected the mathematical aspect the problem, and, general, found good 
agreement between the experimental and analytical results. The earlier works 
Hermann and Keiichi Hayashi’ are notable. method de- 
veloped Krilov® rather convenient and has been introduced into 
book Mr. discusses the same problem most 
complete and thorough manner. this work, addition tables and com- 
plete solutions many special cases, some series solutions are used which 
greatly facilitate the utility the method for practical problems. 

However, the foregoing procedures are thought some difficult 
they involve lengthy mathematical expressions. overcome this difficulty 
method Zusse ASCE, was developed which the entire pro- 
cedure mainly reduced the solution simultaneous equations. The deter- 
mination coefficients and constants many cases may routinized. This 
method involves some legitimate approximations from the engineering point 
view. becomes quite cumbersome, however, the complexity the prob- 
lem increases. 

avoid the tediousness mathematical solution, ingenious electrical 
arrangements have been developed for solving the necessary 
differential equations. The method appears rapid, can provide for the 
variable moment inertia the beam, but requires special equipment and 
techniques. 

overcome some the difficulties mentioned this brief review the 
available methods approximate procedure developed this paper. The 
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method proposed utilizes the great power graphical integration. many 
practical problems the moment inertia beam well the foundation 
modulus are variable quantities along the length the beam. Such variation 
may discontinuous. Even the original mathematical formulation for the 
moment inertia and the foundation modulus may require some approxima- 
tions. Thereafter the solution may become very complex entirely untract- 
able. Although lacks extreme accuracy, the method proposed overcomes 
these complexities. Essentially the problems may solved with nearly the 
same ease whether they are simple complex the sense mathematical 
formulation. The steps used the solution the problem are all very simple 
and are well known engineers. 


GENERAL PROCEDURE 


The usual assumption rigid beam elastic foundation shown 
Fig. 1(a), the true situation for flexible beam shown Fig. and 
idealized elastic foundation only the directly under the beam are 


Original Level 


(b) Elastic Beam Elastic Foundation. 


(c) Elastic Beam Idealized Elastic Foundation. 


Original Level 


(d) Exaggerated Deflections. 


loaded and deflect shown Fig. 1(c). The beam, being flexible, assumes 
some such curve shown exaggerated Since the foundation 
elastic, point the force exerted the spring the beam which 


A 
(a) Rigid Beam Elastic Foundation. 
Original 
Original Level 
A 
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the foundation modulus point and the deflection the beam 
Thus, knowing the foundation moduli, the soil pressure diagram can easily 
constructed manner corresponding deflected beam. Statics considera- 
tions require that the resultant the foundation pressure equal and opposite 
the forces applied the beam, including its own weight. this method 
this requirement must kept mind constantly used over and over 
again adjust the solutions obtained. 

Initially the shape the elastic beam not known and must first as- 
sumed perfectly rigid. This the very assumption used ordinary design 
and its accuracy for the final solution has been questioned herein. Under this 
assumption symmetrically loaded beam will settle equal amount every- 
where the foundation constant uniform variable, and symmetrical 
relation the beam. the loading, the foundation modulus, unsym- 
metrical, the beam will rotate well settle. The amount this 
angular movement can computed easily. Thus (thinking for moment 
terms stresses, which can converted into deflections once the foundation 
modulus known), when there rotation, the stresses ares P/A Mc/I 
instead P/A the former case. The symbols these relations have 
the usual meaning: the resultant the applied loads; simply the area 
under the footing uniform, the transformed area (see Example sub- 
sequently) variable; the distance from the neutral axis the part 
the foundation investigated; the moment inertia the footing contact 
area (transformed area variable) about its neutral axis; and the 
moment, which equal multiplied the distance from the 
centroid the footing area. this manner, regardless the type loading 
the foundation moduli, the pressure distribution under rigid beam may 
and, knowing the pressure distribution, the corresponding deflec- 
tions can determined desired. 

The result such calculations could represented straight line 
Fig. 1(d), and the corresponding pressure distribution under the same 
beam would yz. The subscript this relation designates the particular 
modulus and the particular deflection distance from the end the beam. 
Then the beam considered flexible and its elastic curve determined under 
the influence the actual applied loads and the This done 


plotting the usual subdividing into convenient number 


segments, and graphically arithmetically determining the elastic curve the 
conjugate beam method. (In the expression M/(E I), the symbol denotes 
bending moment section, and the moment inertia the corresponding 
cross-sectional area.) These deflections variously increase and relieve the 
foundation pressures, upsetting the static balance. the beam 
the curved shape just determined the results are adjusted bring them into 
static balance—that is, the resultant the foundation pressure under the 
deflected beam made the equilibrant the applied forces. This pro- 
cedure similar the original treatment given rigid beam, except that 
constant force caused the deflection the beam enters into the calculations. 
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ELASTIC FOUNDATIONS 
Such results are again shown diagrammatically Fig. 1(d) the curved line 
The shape the curve the elastic curve previously deter- 
mined the conjugate beam method. 

Curve exaggerates the deflection the elastic curve was based 
the loading that would occur infinitely rigid beam. the case illus- 
trated, the ends the beam actually would such heavy loading 
deflecting result the flexibility the beam. These deflections, and 
hence the corresponding pressures, the true values. The correct 
answer lies somewhere between the line and the curve Usually 
the best estimate that can made the end this cycle that the true 
elastic curve the average between A’B’ and For the more rigid mem- 
bers the answer closer but for flexible members, closer A’B’. 
the latter case, taken the approximation for subsequent work, 
the solution may not converge. Thus for uniformity and avoid the danger 
working with nonconvergent solution, that the average 
A’B’ and used. 

After averaging the ordinates curves A’B’ and the first cycle the 
solution complete. The corresponding foundation pressures are 
which the average ordinate the two curves point order im- 
prove the accuracy the solution, the pressure distribution corresponding 
the foregoing average deflection then used. This pressure distribution, to- 


gether with the actual loads, used obtain new -diagram, and, pro- 


ceeding before, new elastic curve obtained. Then, adjusted for 
static compatibility. beam “locked” into the curved shape the new 
elastic curve used such calculations. The average the initial elastic 
curve this cycle with the new adjusted elastic curve the desired answer for 
the two cycles. 

has been found that solution made the foregoing basis converges 
very rapidly. Seldom are more than two cycles required. 

illustrate the details this method three examples follow, which some 
pertinent points are discussed. Cases variable moment inertia, variable 
foundation moduli, and the utility the transformed areas for given case 
are shown. facilitate comparison with other methods, data for Examples 
and are from published works cited. However, the type problem 
Example which unusually simple solve the method discussed, nor- 
mally avoided, presents some mathematical difficulties. 


EXAMPLE 


Consider continuous concrete footing having cross section (Fig. 2(a)) 
with line load 150 kips per (k/ft). The elastic modulus concrete 
taken 432,000 kips per ft, and the foundation modulus assumed 
300 kips per ft. this example the soil pressure distribution 
found for the foregoing conditions. The weight the beam neglected. 

Assuming the footing rigid, the pressure distribution 150/30 
kips per and the uniform deflection 5/300 0.0167 ft. 
Considering strip the footing unit width, the bending moment diagram 
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computed for values kips per up, and 150 kips per down (Fig. 
Fig. 2(c), plot the moment inertia the cross-sectional area shown; 


this was computed 5-ft intervals using the relation For ex- 
ample, the center 5.33 ft* per ft, and directly below 


this, Fig. 2(d), the -diagram shown. Only four values are computed 


and the actual curve approximated series straight lines. Only one 
half the diagram shown this problem deals with complete symmetry. 
The M/I-diagram gives the elastic loading for the conjugate beam. The dis- 
tributed loadings are concentrated their respective centers gravity. 
will recalled that such centers for triangles are located distance one 
third from the base. From trapezoid, one way finding the center gravity 
bisect distances such and (Fig. 2(d)) and also lay them off 
the opposite sides shown. This procedure yields points and 
The intersection lines and gives the desired point. 

compute the elastic weights the conjugate beam (Fig. 2(e)) simply add 


the areas the -diagram for particular segment. For example, 


(27.7 70.0) 244 kips per ft. Then the elastic weights are plotted 
vertically scale Fig. 2(f), pole selected, and lines parallel the 
resulting rays are drawn Fig. 2(g). The funicular polygon (Fig. 2(g)), 
certain scale, represents the deflected beam subjected the assumed loading. 
order obtain actual deflections the vertical distances inches are multi- 
plied scale factor. This scale factor obtained multiplying the scale 
the horizontal beam distances the scale the vertical measurements 
the force polygon, the pole distance inches, and, the present example, 
also dividing the elastic modulus the beam. Thus, the factor the 
ordinates from straight line passing through the end points the beam are 
recorded Fig. 2(g). The pole distance may adjusted make the scale 
factor more convenient. 

The foundation modulus 300 kips per The deflections, feet 
(Fig. 2(g)), correspond the pressures (Fig. 2(h)) which are obtained the 
relation provided that, for awhile, the ends are assumed precisely 
the surface. Then the beam with the curvature shown and the 
entire curve depressed uniformly into the foundation obtain static bal- 
ance. Summing the foundation pressures already developed when the ends 
kips. Hence, 150 kips remains uniformly distributed; 


original solution was 


that 2.03 kips per ft. Adding this value the ordinates 


Fig. 2(h) the results are shown Fig. Averaging this with the initially 
assumed uniform pressure distribution kips per yields the results 
one cycle the semi-graphical solution. This shown Fig. 
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improve the accuracy the results, another cycle may made. Soil 
pressure loading shown Fig. together with the actual loads, are used 


the next approximation. The -diagram shown Fig. 2(k), and the 


elastic curve, Fig. the accompanying force polygon not shown. After 


TABLE 1.—REFINEMENT Pos- 
SIBLE 
EXAMPLE (PRESSURE 
INTENSITY Kips 
PER SQUARE 


adjusting the results for static compati- 
bility that accommodate the curve (Fig. 
the results shown Fig. 2(m) 
are obtained. Averaging these results 
with those Fig. gives the results 
the second cycle (Fig. 2(n)). 


Two cycles are believed give re- 
Rigid First Third that are sufficiently accurate for 
cycle cycle cycle 


segments must used, and, perhaps, 


few more cycles work must per- 

formed. comparison the results, 
the foregoing two cycles and one additional cycle are shown Fig. 2(0) and 
Table 


EXAMPLE 


“weightless” beam, in. in., with the loading shown Fig. 3(a) 
resting elastic foundation. The modulus the foundation 200 
percuin. The elastic modulus the beam 1.5 per in. (psi), 
and its moment inertia Pressure distribution under this beam 
determined this point. (Mr. Hetényi has solved this problem mathema- 
Data for Example were selected from his work.) 

The resultant the applied loads found located 7.5 in. off center 
shown Fig. This indicates that, addition the uniform settlement, 
angular displacement takes place. The contact area the beam with the 
foundation 1,200 in. and the moment inertia this area 
1,440,000 Hence, concluded that, the beam were rigid, 
the pressure end would —11.23 per in., and end —5.11 per 

3.06. This pressure variation linear shown Fig. 3(c). usual, 
minus signs indicate compression. 

With actual applied loads Fig. and pressure distribution 
Fig. 3(c), the bending moment diagram Itis approximated with 
the series straight lines shown Fig. 3(d). After computing the elastic 
weights, plotting them suitable scale, and selecting pole (just the 
previous example), the elastic line shown Fig. 3(e) obtained. The de- 
flections from the line passing through the ends are noted. Lack symmetry 


this case required that the entire curve drawn in. The scale factor had 


Elastic Foundation,’’ Hetényi, The Univ. Michigan Press, Ann Arbor, Mich., 
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divided the moment inertia the beam, bending moment rather 
than the -diagram was used. 


the beam the deflected shape just found, ends and are 
moved that they are the surface. The pressures developed the founda- 
tion the curved beam this position These are shown Fig. 3(f). 
The resultant these pressures 12.15 kips (12.15*) upward, 1.8 in. the 
left the center line. These quantities were obtained subdividing the 
pressure area into triangles and rectangles and taking their moments about 
one end the beam. This resultant excess the applied loads which 
total 918 kips; therefore, necessary apply upward pull 2.35 kips 
the curved beam. Moreover, the pressure resultant develops 
clockwise resisting moment only 21.8 kip-in., whereas, keep the applied 
forces balance, 73.5 kip-in. moment required. develop the founda- 
tion resistance for the difference 51.7 kip-in., additional rotation the 
curved beam (as unit), this magnitude, must applied the counter- 
clockwise direction. Hence, from the known relations these pressure increments 


in.; and, similarly, persq in. This pressure distribution varies 


linearly from point point 


TABLE Adding algebraically the pres- 


(PRESSURE INTENSITY sure distribution shown Fig. 3(f), 
PER SQUARE the distribution shown Fig. 3(g) 
obtained. This diagram stati- 


cally compatible with the applied 


9.70 109 11.2 the applied loads. 

Center* 10.5 10.7 10.39 The average the pressure 


Pressures the center are obtained inter- the desired results for one 
polation. shown Fig. 3(h). 
For the second cycle, the pressure 
distribution Fig. 3(h), together with the applied loads, gives approximate 
moment diagram and the new elastic curve Figs. and respectively. 
this curved beam and adjusting the results for static compati- 
bility, the pressure distribution shown Fig. obtained. The average 
the values Figs. 3(h) and 3(k) represents two cycles (see Fig. 

For comparison, the results one additional cycle and the results 
mathematical solution,“ are given Table and Fig. 3(m). The agreement 
appears satisfactory. 


EXAMPLE 


Consider weightless” beam, in. in. (Fig. 4(a)), with two equal loads 
100 each the ends. The constants for the beam are: 2.5 
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per in. and 144 The foundation modulus assumed vary 
end Mathematical solutions for the deflection line and pressure distribution 
for these data are Herein solved the semi-graphical method. 


0.90 Pounds per inch 


45 


per inch 


3 


58 
2.04 
Pounds per inch 


‘ 2 3 4 s é 7 


This problem differs from the previous cases the foundation modulus 
varies along the length the beam. The stiffer foundation end offers 
more resistance the load. However, precisely the same amount resistance 
would obtained with uniform foundation modulus but with the width 
the beam increased decreased the ratio the foundation modulii. the 
example considered, designated the uniform base modulus, the width 
the beam end becomes in. the foundation modulus, 
this case, varies linearly, the equivalent area over the region where the 
modulus trapezoid shown Fig. 4(b). course, any other modulus 
may used base. The transformed area the contact surface 3,840 
in.; its centroid in. from end and the moment inertia this 
area about its centroidal axis 3,744,000 

The loading symmetrical this case but the center the foundation 
pressure resistance not. The remaining parts the problem are similar 


Elastic Hetényi, The Univ. Michigan Press, Ann Arbor, Mich., 
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Example The pressure distribution for rigid beam point is: 


200 


3,744,000 
per in. the actual rigid beam, in. being the width the trans- 
Extra ner pressure distribution therigid beam 


EXAMPLE (PRESSURE has been computed 20-in. intervals. 
INTENSITY Results are plotted Fig. 4(c). This 
PER pressure distribution, together with the 


two loads, yields the moment diagram 


1.21 0.95 dation, the pressure distribution caused 
0.90 


179 The deflections Fig. 4(f) are 
multiplied the corresponding variable 
foundation modulii. The resultant this pressure, together with the two end 
loads, gives total force 602.4 acting 9.8 in. the right the centroidal 
axis. This adjusted statically compatible, using the transformed con- 
0.275 2.20 per in. actual beam, etc. These results are shown 
Fig. 4(g). Averaging this curve with the initial distribution shown Fig. 4(c), 
the results for one cycle are shown Fig. 4(h). Results for two cycles 
together with the analytical are given Table and Fig. Again, 
the agreement found may considered entirely satisfactory. 


tact area: 


The paper, with the illustrative examples, treats most the important 
points that are encountered practice connection with beams elastic 
foundation. Rather complex problems with variable moments inertia 
the beam, and variable foundation moduli can solved with relative ease. 

the examples, the writer assumed consistently that the foundation can 
develop tension well compression stresses. This detail importance 
during the performance calculations. However, the end results must 
interpreted with caution. soil problem, tension possible. The 
part the beam that experiencing only upward pressure distribution must 
statically compatible. Likewise, must borne mind that, since ex- 
amples involving soils were used, the foundation modulus not simple 
quantity. Reaction does not depend only deflection point; the ad- 
joining areas also affect this value. provide for this influence adjoining 
areas, however, one the merits the suggested procedure may com- 
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pensated the light this knowledge) and the problem remains comparatively 
simple. 

With minor modifications, the method discussed well adapted other 
similar problems. Its utility for flexible floating structures particularly 
notable. 

The development the paper followed the semi-graphical sequence; but the 
method similar the well known moment-area procedure, and, desired, 
those who abhor graphics, may set tables and perform all the necessary 
calculations accordingly. attempt was made this paper illustrate the 
power and simplicity the graphical method. All procedures are conventional. 
should require little time master the over-all scheme. 

The method not proposed for use situations which ready analytical 
solutions are available. Unfortunately, these cover only the simple cases. 
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